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The operation of delumping is based on the fact that the oil reservoir has only two non aqueous phases 
(liquid and gas) , allowing fluid regrouping and detailed decomposition calculations to be carried out. 
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METHODE DE SIMULATION POUR PREDIRE EN FONCTION DU TEMPS UNE COMPOSITION DETAILLEE 
D'UN FLUIDE PRODUIT PAR UN RESERVOIR. 



_ Methode de simulation pour predire en fonction du 
temps une composition detaillee a Q composants et/ ou 
pseudo-composants d'un fluide produit par un reservoir. 
Elle est du type consistant a: 

- representer le reservoir sous la fornne d'un reseau de 
mailles 0) dont chacune constitue un volume elementaire 
rempll de fluide, 

- realiser, de fagon connue en soi, une simulation com- 
positionnelle du fluiae regroupe a N composants et pseudo- 
composants (i), N etant inferieur a 0, et elle est caracterisee 
en ce qu'elle consiste en outre a: 

e) determiner, a chaque pas le temps (m) et pourcha- 
que maille (j) la composition des phases liquide et vapeur du 
fluide detaille de la composition 0 du fluide, a partir des va- 
leurs de la fraction vaporisee et des constantes d'equi- 
libre (kj^ du fluide regroupe, 

f) evaluer pour chaque maille, au pas de temps (m+1), la 
quantite molaire de chacun des Q composants et/ ou pseu- 
do-composants du fluide detaille a partir des valeurs con-es- 
pondantes au pas de temps (m), des debits des phases du 
fluide regroupe et des compositions des phases liquide et 
vapeur du fluide detaille determlnees a Tetape e). et a 

g) evaluer, pour chaque puits de production, la compo- 



sition detaillee du fluide produit entre des instants t et t' cor- 
respondant aux pas de temps m et m+1 , a partir des debits 
des phases du fluide regroupe et des compositions des pha- 
ses liquide et vapeur du fluide detaille determlnees a I'etape 
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Methode de simulation pour predire en fonction du temps 
une composition detaillee d'un fluide produit par un reservoir 



La presente invention conceme nne methode de simulation pour 
predire, en fonction du temps, la composition detaillee d'un fluide produit 
par un reservoir et plus particulierement la composition detaillee d'un fluide 
contenu dans et produit par un gisement petrolier dans lequel sont implantes 
im ou plusieurs puits de production. 

La simulation compositionnelle d*un gisement petrolier est 
couramment utilisee pour foumir les profils previsionnels de production du 
gisement qui permettent notamment de determiner le schema de production 
le mieux adapte a ce gisement 

La simidation compositionnelle d'un gisement petrolier est mise 
en oeuvre non pas en utilisant une description reelle du fluide mais en 
utilisant un fluide modelise par un nombre de composants plus reduit que le 
nombre de composants du fluide reel. En effet, le nombre de composants du 
fluide de gisement etant relativemeut grand, la modelisation avec tons les 
composants conduirait a des temps de calcul trop importants. Ce temps de 
calcul prohibitif a amene les specialistes a regrouper les composants purs en 
pseudo-composants, par exemple des pseudo-composants regroupant Tazote 
et le methane, des pseudo-composants regroupant les hydrocarbines en C3 
et C4, ... etc., et a realiser la simulation de gisement sur un nombre reduit 
N de composants purs et de pseudo-composants. Une telle composition 
reduite est designee par Texpression "composition regroupee du fluide" ou 
"flxiide regroupe". Generalement N represente 5 a 8 composants et pseudo- 
composants, ce qui est considere comme etant sufFisant pour bien 
representer le comportement du fluide de gisement aux conditions de fond. 

Avant d'effectuer la simulation compositionnelle de gisement, on 
represente le reservoir ou gisement sous la forme d*un reseau de mailles 
dont chacune constitue un volume elementaire dudit reservoir. Le nombre 
de mailles peut atteindre plusieurs milliers et chaque maille presente des 
proprietes qui lui sont propres, comme la geometric, la porosite, la 
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permeabilite, etc. . . De plus, a certaines mailles peut correspondre au moins 
un puits d'injection ou de production qui est implante dans ledit reservoir. 

La simulation de gjsement permet de calculer pour chaque maille 
un certain nombre de variables principales du fluide de gisement, lesdites 
5 variables pouvant etre la quantite (nombre de moles) de fluide, la compo- 
sition ou fraction molaire de chaque composant et pseudo-composant du 
fluide regroupe et la pression regnant dans chaque maille. Ces variables 
sont connues a Tinstant t=0 (debut de Texploitation du gisement) et sont 
ensuite calculees par la simulation de gisement a chaque instant t ou pas de 

10 temps m. Pour chaque maille, il est possible, a partir de ces variables 
principales, de calculer toute autre propriete du fluide present dans la 
maille, comme le nombre de phases, la composition de chaque phase, etc. 
Pour chaque maille a laquelle est associe un puits de production, il est 
possible, a partir de ces variables principales et des contraintes de 

15 production imposees, de calculer egalement le debit de production et la 
composition du fluide prodxiit par ledit puits. 

Les proprietes thermodynamiques du fluide peuvent etre 
calculees en utihsant une des equations d'etat bien connues des specialistes 
et qui ne seront done pas rappelees. 

20 La modelisation a N composants, si elle est suffisante pour 

representer le comportement du fluide de gisement aux conditions de fond, 
n'est plus appropriee poxu" simuler le comportement du fluide dans les 
installations pour son exploitation en surface, laquelle exploitation necessite 
de connaitre une composition plus detaillee a Q composants et/ou pseudo- 

25 composants du fluide produit par le gisement, Q etant superieur a N et par 
exemple de I'ordre de 20 a 30. 

Jusqu'ici, on obtenait cette composition detaillee a Q composants 
et/ou pseudo-composants a partir de la composition regroupee N en 
supposant que la composition des pseudo-composants restait constante dans 

30 le temps. 

Une telle maniere de proceder est source d'erreurs car la 
composition des pseudo-composants varie dans le temps, notamment en 
fonction de la pression du gisement ou suite a Tinjection d'un gaz dans le 
gisement. 
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De ce fait, la composition detaillee obtenue reste tres approxima- 
tive et ne permet pas de predire de maniere satisfaisante le comportement du 
fluide dans les conditions d'exploitation en surface bien des annees apres le 
debut de la production. 
5 La presente invention a pour but de proposer une methode de 

simulation pour predire la composition detaillee du fluide produit par im 
reservoir en tenant compte des parametres du reservoir et qui soit beaucoup 
plus precise. 

L'invention a pour objet une methode du type consistant a : 
10 a) representer le reservoir sous la forme d*un reseau de mailles (j) 

dont chacune constitue xm volume elementaire rempli de fluide, 

b) definir le fluide par une modelisation regroupee a N 
composants et pseudo-composants (i) et a determiner xme equation d'etat 
decrivant le fluide dans cette modelisation regroupee, 
15 c) definir egalement le fluide par une modelisation detaillee a Q 

composants et/ou pseudo-composants, Q etant superieur a N, et a deter- 
miner une autre equation d'etat decrivant le fluide dans cette modelisation 
detaillee, 

d) realiser, de fa9on connue en soi, une simulation composi- 
20 tionnelle du fluide regroupe a N composants et pseudo-composants (i), 

ladite simulation compositionnelle pennettant de calculer au moins pour 
chaque maille (j) et a des pas de temps consecutifs (m, m+l, ...) la fraction 
vaporisee (e^), les constantes d'equilibre liquide-vapeur [k"^) de chaque 

composant (i), les debits d'injection ou de production (s") et pom* chaque 
25 paire de mailles (j, h) les debits des phases liquide (u^,^) et vapeur (u^j,) du 

fluide a N composants et pseudo-composants, et elle est caracterisee en ce 
qu*elle consiste en outre a: 

e) deteraiiner, a chaque pas de temps (m) et pom chaque maille 
(j) la composition des phases hquide et vapeur du fluide detaille de la 

30 composition Q du fluide, a partir des valeurs de la fraction vaporisee (9^) et 

des constantes d'equihbre {k^) du fluide regroupe, 

f) e valuer pour chaque maille, au pas de temps (m+1), la quantite 
molaire de chacun des Q composants et/ou pseudo-composants du fluide 
detaille a partir des valeurs correspondantes au pas de temps (m), des debits 
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des phases du fluide regroupe et des compositions des phases liquide et 
vapeur du fluide detaille determinees a I'etape e), et a 

g) evaluer, pour chaque puits de production, la composition 
detaillee du flxiide produit entre des instants t et t" correspondant aux pas de 
5 temps m et m+l, a paitir des debits des phases du fluide regroupe et des 
compositions des phases liquide et vapeur du fluide detaille determinees a 
Tetape e). 

Selon une autre caracteristique de la presente invention, les 
etapes e, f et g sent raises en oeuvre en meme temps que la realisation de la 
10 simulation compositionnelle du flxiide regroupe. 

Selon une autre caracteristique de la presente invention, les 
resultats de Tetape d sont memorises puis utilises ulterieiu-ement pour la 
mise en oeuvre des etapes e, f et g. 

La methode selon la presente invention est mise en oeuvre pom- 
IB un reservoir constitue, par exemple, par un gisement petrolier. De maniere 
connue, on realise xme simulation de gisement. A cet effet, on represente le 
gisement sous la forme d*im reseau maiUe, certaines des mailles ou un 
groupe de mailles etant associes a un des puits de production rmplantes dans 
le gisement petrolier a exploiter. 
20 Pour reduire le temps de calcul, la simulation compositionnelle 

du gisement est effectuee sur un nombre limite N de composants et pseudo- 
composants, par exemple de 5 a 8, definis de la maiuere indiquee 
precedemment, ces composants et pseudo-composants etant selectionnes en 
fonction de la nature du gisement. 
25 Le fluide regroupe a N composants et pseudo-composants est 

decrit par xme equation d'etat qui pourrait etre par exemple celle de PENG- 
ROBINSON adaptee a ladite composition regroupee. 

Au debut de Texploitation du gisement, on cormait, par des 
mesures effectuees au prealable, pour chaque maille j du reseau, les 
30 variables principales a Finstant t==0 ou pas de temps m=0. Ces variables 
principales sont : 
la quantite de fluide : f 

la fraction molaire du composant i : z" avec 1 < i < N 
la pression : p^ 
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Pour passer du pas de temps m au pas de temps m+1, la 
simulation compositionnelle comprend plusieurs etapes detaillees ci-apres. 

Dans une premiere etape et pour le pas de temps m, un calcul 
d'equilibre liquide-vapeur (flash) est effectue sur le fluide regroupe, pour 
5 chaque maille j afin de determiner : 

la fraction vaporisee 9 J" 

la fraction molaire de chaque composant i dans la phase huile (x" avec i = 1...n) et 
la fraction molaire dudit composant i dans la phase gaz (y^ avec i = 1. ..n), 

Ces diverses fractions sont determinees avec le systeme 
10 d'equations suivant : 

<I>.(x,7 X- ) et 4>,(y7,..-.yK,) avec 1 < i < N 

et 

N N 

1=1 1=1 

15 

dans lesquelles : 

Oi est la ftigacite du composant i. 

Selon une premiere caracteristique de I'invention, on memorise la 

fraction vaporisee et les constantes d'equilibre k"^ - avec i variant 

20 entre 1 et N, afin de les utiliser dans une autre etape denommee 
"delumping" qui sera explicitee ulterieurement. 

Dans une deuxieme etape et pour chaque maille, on evalue 
certaines des proprietes de chaque phase liquide et gaz, a savoir masse 
volumique (p^ et p^), viscosite, saturation et permeabilite relative. 

25 Dans une troisieme etape, pour chaque maille j et au pas de 

temps m, on calcule les coefficients de Tequation poxir la pression. 

Les equations pour la pression d'une maille j a I'instant t+ At ou au 
pas de temps m+1, sont : 



2775094 



I 



20 



^ At J heJ(j) 



Pour resoudre ces equations, les coefficients a calculer sout 
la compressibilite totale Cj^ (fluide + roche), 

les transmissivites generalisees 

Trh = Z<(^oXi+\yJf-t-, 



10 - le terme source • sj" = ^ 

1=1 

Dans ces equations : 
h sont les mailles voisines de la maille j, h e J(j), 

a™ est le volume molaire partiel du composant i dans la maille j au pas de 
temps m, et X,^ sont les mobilites des phases huile (o) et gaz (g), 
15 tjt, est la transmissivite entre les centres des mailles j et h, 1 etant la maille 

amont. 

Dans xme quatrieme etape, on resout pour toutes les mailles les 
equations pour la pression a Tinstant t+At c'est-a-dire au pas de temps m+1. 
Ces equations forment un systeme lineaire du type ; 



A. p = B 



ou p = (pI"*^ ,...pP*)^, J etant le nombre de mailles du reseau. 

Ces equations sont resolues par un procede iteratif standard. 
25 Dans une cinquieme etape, on calcule au pas de temps m+1, les 

debits des phases entre chaque paire de mailles adjacentes et pour chaque 
puits associe a partir, par exemple, de la loi de DARC Y. 

u:h=-p:?^oTtjH(pr'-pr') 

30 
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1 etant la maille amont avec l=j ou h. 

Selon une autre caracteristique de rinvention, on stocke ou 
memorise les debits des phases qui seront utilises dans le delumping. 

Dans une sixieme etape, on evalue la quantite f et la composition 
5 z du fluide regroupe, dans chaque maille j et au pas de temps m-f 1 a partir 
des valeurs au pas de temps m, des debits des phases evalues a la cinquieme 
etape et des compositions des phases calculees a Tetape 1 . Pour ce faire, on 
applique les lois de conservation definie par les equations suivantes : 

f;*' = f;-At|2te+u-J-sj 

10 Par un processus iteratif, on repete les etapes ci-dessus pour 

calculer les variables principales aux temps t+2At, t+SAt, etc... ou encore 
aux pas de temps m+1, m+2, m+3, etc., I'increment At pouvant etre 
constant ou variable. 

Selon la presente invention, on definit une modelisation detaiUee 

15 a Q composants et/ou pseudo-composants ainsi qu'une autre equation d'etat 
decrivant le fluide ainsi modelise. 

Pour cela, on choisit prealablement les composants et/ou pseudo- 
composants dont on veut connaitre les caracteristiques et qui sont 
necessaires pour representer le fluide de production dans les conditions de 

20 surface. Le nombre Q qui est superieur au nombre N de la composition 
regroupee est generalement de Fordre de 20 a 30. On connait la composition 
du fluide detaille a Q composants et/ou pseudo-composants au pas de temps 
m=0. 

Les donnees stockees aux premiere et cinquieme etapes sont 
25 utilisees pour les etendre aux Q composants et/ou pseudo-composants de la 
composition detaiUee. 

Conune Tecoulement pent s'effectuer dans n'importe quelle 
direction entre les mailles voisines, chaque maille ne peut pas etre traitee 
independamraent, C'est pourquoi, a chaque pas de temps, toutes les mailles 
3 0 sont traitees ensemble. 
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Pour reaiiser le delumping, on fait I'hypothese qu'a chaque pas de 

temps m, 

- la jfraction molaire de la phase vaporisee dans chaque maille est 
independante du nombre de constituants du fluide. 



0^ = 97 



- le nombre de moles de liquide et de vapeur s*ecoulant entre les mailles et 
dans ou hors des puits est independant du nombre de constituants du 
fluide. 



10 



Les lettres en majuscule sont utilisees pour representor le fluide 
detaille et les lettres en minuscule sont utilisees pour representer le fluide 
regroupe. 

Le calcul des constantes d'equilibre peut etre effectue a partir de 
15 Texpression : 

Ln[K,]- ACo -hX^C^.Ti. ou AC„ = -C^ 

comme indique par C. LEIBOVICI dans son article "A Consistent 
Procedure for Pseudo-Component Delumping" pam dans FLUID PHASE 
20 EQUILIBRIA, 117 (1996), 225-232. 

Pour une equation d*etat a deux parametres, comme ceUe de 
PENG ROBINSON, I'equation generalisee de LEIBOVICI peut etre ecrite 
' pour le fluide detaille de la maniere suivante : 

25 Ln(K,7) = <+P7V^" + Yrbi 

ou ai et bi sont les parametres de Tequatioa d*etat pour le composant I et ou 
les parametres ti^,P^ ety" sont calcules a partir des constantes d'equilibre 

pour le fluide regroupe {k^) en minimisant la fonction : 

1=1 ^ 
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De ce fait, lorsqu'on connait la composition detaillee Zi7 dans la 
maille j au pas de temps m, on peut estimer les compositions detaillees des 
phases liquide et vapeur dans la meme maille j et pour le meme pas de 
temps m. 



7"* K'^Z™ 
X;"- 7^ r Yr = 7-^^ r avecl<I<Z 

i+07(k;^-i) ' i+0;"(k^-i) 

I etant le composant et/ou pseudo-composant de la composition 
detaillee. Ces compositions peuvent etre nomialisees afm que leur somme 
10 soit egale a 1 si necessaire. 

Lorsqu'on coimait les compositions des phases detaillees pour 
toutes les mailles du reseau au pas de temps m, il est possible destimer les 
compositions detaillees ziT' au pas de temps consecutif m+1 a partir des 

equations ci-apres : 



ZS*Fr-(Yl?So^ + Xf?S™)- E (Yi?U^h + Xreh) 
Zr' = avecI=l,...,Q 



ou 



20 Fr'=Fr-(s,";+s^)- Eta+uSh) 

dans lesquelles j -j pour un debit positif de la maille j vers la maille h ou 
dans le puits et j-h pour un debit negatif, j' correspondant au fluide injecte 
pour les puits d'injection on S sera negatif. 

Grace a la presente invention, on obtient une composition 

25 detaillee a Q composants et/ou pseudo-composants qui permet de nueux 
definir les profils de production previsionnels et ce, avec une precision 
accrue comparee aux methodes anterieures. Cette plus grande precision des 
profils de production previsionnels permet un choix de schema de 
developpement (type et taille dinstallations de surface, nombre de puits, 

30 etc..) plus liable, avec pour consequence xme augmentation de la rentabilite 
et une diminution du risque economique. 
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Dans ce qui precede, certaines valeurs obtenues au cours de la 
simulation de gisement sont stockees ou memorisees pour etre utilisees 
ulteriexxrement dans le delumping. II va de soi que ces memes valeurs 
pourraient etre utilisees directement si on decide d'effectuer le delumping en 
5 meme temps que ladite simulation de gisement. 
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REVENDICATIONS 



5 1. Methode de simulation pour predire, en fonction du temps, la 

composition detaillee d'un fluide produit par un reservoir, du type consistant 
a : 

a) representer le reservoir sous la forme d'un reseau de mailles (j) 
dent chacune constitue un volume elementaire rempli de fluide, 
10 b) defmir le fluide par une modelisation regroupee a N compo- 

sants et pseudo-composants (i) et a determiner une equation d'etat decrivant 
le fluide dans cette modelisation regroupee, 

c) definir egalement le fluide par une modelisation detaillee a Q 
composants et^ou pseudo-composants, Q etant superieur a N, et a deter- 

15 miner une autre equation d'etat decrivant le fluide dans cette modelisation 
detaillee, 

d) realiser, de fa<;on connue en soi, une simulation composi- 
tionnelle du fluide regroupe a N composants et pseudo-composants (i), 
ladite simulation compositionnelle permettant de calculer au moins pour 

20 chaque maille (j) et a des pas de temps consecutifs (m, m+1, ...) la fraction 
vaporise (ej*), les constantes d'equilibre liquide-vapeur [k^) de chaque 
composant (i), les debits d'injection ou de production (s^) et pour chaque 
paire de mailles (j, h) les debits des phases liquide (u,";J et vapeur (u^J du 
fluide a N composants et pseudo-composants, et elle est caracterisee en ce 

25 qu'elle consiste en outre a : 

e) detenmner, a chaque pas le temps (m) et pour chaque maille 
(j) la composition des phases liquide et vapeur du fluide detaille de la 
composition Q du fluide, a pardr des valeurs de la fraction vaporisee 

des constantes d'equilibre [k"^) du fluide regroupe, 
30 f) evaluer pour chaque maille, au pas de temps (m+l), la quantite 

molaire de chacxm des Q composants et/ou pseudo-composants du fluide 
detaille a partir des valeurs correspondantes au pas de temps (ra), des debits 
des phases du fluide regroupe et des compositions des phases liquide et 
vapeur du fluide detaille determiaees a I'etape e), et a 
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g) evaluer, pour chaque puits de production, la composition 
detaillee du fluide produit entre des instants t et f coirespondant aux pas de 
temps m et m+1, a paitir des debits des phases du fluide regroupe et des 
compositions des phases liquide et vapeur du fluide detaille determinees a 
5 I'etape e). 

2. Methode selon la revendication 1, caracterisee en ce que ies 
etapes e, f et g sont mises en oeuvre en meme temps que la realisation de la 
simulation compositionnelle du fluide regroupe, 

3. Methode selon la revendication 1, caracterisee en ce que les 
10 resultats de Tetape d sont memorises puis utilises ulterieurement pour la 

mise en oeuvre des etapes e, f et g. 
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